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Abstract 

The vortex state in a single crystal of the layered organic superconductor k-(BEDT- 
TTF2Cu(NCS)2, where BEDT-TTF (or FT) is bis(ethylenedithio)tetrathiafulvalene, was studied 
by ^H-NMR. Under a low field region around 0.75 T, the vortex glass-liquid transition was demon- 
strated by a diverging of the longitudinal nuclear spin relaxation rate and peak-broadening in 
spectra. Under a high field region near the upper critical field Hc2{0) —7 T, the curvature of 
nuclear spin relaxation curves showed a drastic change at the temperature where the emergence 
of the quantum vortex slush state was reported. The mechanism in this curvature change was 
discussed in terms of the fluctuating field produced by fragments of vortex glass. 
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I. INTRODUCTION 



Quantum fluctuation often brings a non-trivial physical state to atoms or solids at low 
temperatures, represented by liquid ^He or spin-gap systems. The former persists to be liquid 
even at the absolute zero, and the latter shows a non-magnetic spin liquid state, where the 
Neel state is destructed by the quantum spin fluctuation. Quantum melting of the flux lattice 
in superconductors has attracted much interest, because it is expected to bring about novel 
and exotic states of vortices. In superconductors with a quasi-two-dimensional structure, a 
short inter-plane spacing and a high critical temperature, the fluctuation is expected to have 
a large effect on the vortex state when the magnetic fleld is applied perpendicular to the 
conduction plane so that the intrinsic pinning has little effect on the motion of vortices l|, 
a a. Although thermal fluctuation dominates the vortex state in the high temperature 
region near T^, Ikeda[3l] has shown that there is a crossover temperature Tqf, below which 
a quantum fluctuation should be dominant, and that Tqf is proportional to H in the quasi- 
two-dimensional systems, so that the quantum melting phenomenon is expected to exist at 
a high fleld near Hc2- In fact, Sasaki et al. made an intensive study on the layered organic 
superconductor k-(BEDT-TTF)2Cu(NCS)2, which is one of the most anisotropic systems in 
K-type organic complexes, and showed that in the vicinity of iJc2— 7 T, a flnite resistive state 
with non-ohmicity persists at low temperatures down to 100 mK.pJ] They conjectured from 
the results of transport measurements that the observed low resistance region corresponds 
to the slush state, where a short range order of vortices remains. The microscopic nature 
of such slush is still unexplored. The purpose of this study is to investigate the quantum 
vortex slush state microscopically by NMR. 

Figure 1 shows schematic HT-phase diagram of the vortex state in k-(BEDT- 
TTF)2Cu(NCS)2 reported in ref. [ll, where the phase boundary between the vortex glass 
and the vortex liquid, which is usually referred as the irreversibility line H-^-^^. intercepts the 
fleld around 4 T much lower than Hc2- The slush state is located at low temperatures in the 
fleld region higher than 4 T. So far, NMR studies on the title compound in powder samples 
were reported by Takahashi et al. Q, H, Q They found that under a low applied fleld 
around a Tesla, the longitudinal spin relaxation rate shows a pronounced diverging be- 
havior at temperatures far below T^. This anomaly, which coincides with the irreversibility 
line, is now understood to be due to the critical slowing down of magnetic fleld fluctuation 
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FIG. 1: HT-phase diagram of k-(BEDT-TTF)2Cu(NCS)2 reported by Sasaki et al.[l] and the 
experimental conditions of temperature range and fields performed in this paper, denoted by t-wo 
grey thick dotted lines. Each curves sho-w the boundary bet-ween the regions of the normal state, 
the vortex liquid phase, the vortex glass phase, the Bragg glass phase, and the quantum vortex 
slush state. Symbols of data points on curves are the same as in Ref. 



at the vortex glass-liquid transition'15']. Further examination on single crystals under a lo-w 
field -were reported by other groups. Van-quynh et al. observed that the ^H-NMR spectrum 
consists of a narrow peak and a broad peak overlapping each other, and argued that the 
latter comes from the part of the vortex lattice distorted by NMR excitation pulseslTtl. u-SR 

n Q 

studies[8|, |9| report that a ske-wed Abrikosov pattern appears in the muon rotating spectra 
only under an extremely lo^w transverse field (TF) of 25 Oe, corresponding to the Bragg- 
glass phase, and not at 400 Oe, corresponding to the vortex-glass phase. All these previous 
reports dealt with the low field region, where vortices solidify at low temperatures. 

In this paper, we report on a microscopic investigation of the vortex slush state located 
at high field region above 4 T by ^H-NMR, which is a sensitive probe for the magnetic flux 
inside a superconductor. 

II. EXPERIMENTAL 

Single crystals of k-(BEDT-TTF)2Cu(NCS)2 were synthesized with an electrical oxidation 
method. Two crystals jjl, jj2 with an approximate size of 3x1.5x0.1 mm were selected. ^H- 
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FIG. 2: Temperature dependence of Tf ^ and the peak width (FWHM) of ^H-NMR on the sample 
jjl under the low field of 0.75 T, which was applied perpendicular to the conducting 6c-plane. Solid 
and dashed curves are guides for the eyes. The glass melting temperature Tm and the critical 
temperature Tc under Hq is shown by solid arrows. The inset shows typical spectra profiles. 
Horizontal dashed arrows indicate the definition of the line width (FWHM). 

NMR measurements were carried out separately on each two crystal by using a conventional 
^He cryostat set in a 6T cryogen-free superconducting magnet with a field homogeneity of 
50 ppm in 10 mm dsv. In all the NMR measurements, the sample was field-cooled, with a 
slow coohng rate around 10-20 K/hour. Spectra were obtained by superposing each Fourier- 
transformed spectra of the spin-echo beats taken under several magnetic fields with an 
interval of 10 Oe. The spin-lattice relaxation rate Tf ^ in the low field region was determined 
by fitting the dependence of the echo amplitude on the repetition interval r of measurements 
to the function 



where Iq is the nuclear magnetization at thermal equilibrium. In the high field region, where 
the slush state emerges, Ti exceeds a few hundred minutes, making it difficult to determine 
its absolute value. Instead, we have traced the initial curvature of the relaxation curves and 
investigated its temperature dependence. 
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III. RESULTS AND DISCUSSION 



Figure 2 shows the temperature dependence of Tf ^ at a low field of 0.75 T, which 
corresponds to the line (L) in Fig. 1. The field was applied perpendicular to the con- 
ducting 6c-plane. One can see that T^^ diverges around 3 K, which is much lower than 
Tc{H = 0.75T) ~ 8.5 K. This temperature coincides with the vortex liquid-glass transition 



reported by NMR on powder samples 
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14| and by magnetization measurements 



and we can safely assign the divergence in T^^ to the vortex glass melting, where 
the fiuctuation in the magnetic field produced by vortex motion shows a critical slowing 
down. Generally, the nuclear spin-lattice relaxation is driven by a Fourier component of the 
Lamor frequency, that is, around 100-200 MHz, in a fiuctuating magnetic field. This slow 
component is prominently enhanced in the vicinity of second order phase transitions such 
as glass meltingjn, 12, [l4 |. 

The temperature dependence of the width of the spectral peak defined as the full width 
at half maximum (FWHM) is also shown in Fig. 2. With decreasing temperature, the 
width shows a slight increase from 8 to 11 Oe, and an abrupt increase to 14 Oe at the 
temperature slightly lower than that where T^^ diverges. This difference comes from the 
fact that the T^^ diverges when the dominant frequency component of the field produced 
by vortex coincides with the Lamor frequency while the spectrum broadens when the field 
becomes nearly static. 

The observed spectra shown in the inset always consist of a symmetric peak and bear 
neither skewness nor splitting in the entire temperature range, in contrast to the previous 
report 0]. Generally, in the presence of a rigid Abrikosov fiux lattice with translational sym- 
metry, the internal field distribution must show a characteristic skewed pattern, depending 
on the lattice symmetry [l^ . The absence of the skewness in the observed spectra indicates 
that the vortex pattern loses its translational symmetry within each plane. This is consistent 
with the fact that the system resides in the vortex glass state under an applied field of 0.75 
T. 

Next, we proceed to the high field data, the field and temperature region of which cor- 
responds to the line (H) in Fig. 1. The spectrum consists of a single peak also in this field 
region, and its width of 14(5) Oe does not show any temperature dependence between 15 and 
0.5 K. Figure 3 shows the initial part of typical relaxation curves at various temperatures. 
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FIG. 3: Typical relaxation curves at high field around 5T, a) for initial part within 60 s for the 
sample (Jl, and b) within 10^ s of the sample ji2. Solid and dashed curves indicate the power law 
function. 

Under a high field of around 4-5 T, the relaxation time is drastically elongated to be far 
above 10^ s in the measured temperature region below 6 K. Therefore, the absolute value of 
Ti and hence its temperature dependence could not be precisely determined in the allowable 
experiment time. The initial part of the relaxation curves, at least from zero to 2x10^ s, 
exactly follows the power law of time as shown by the log-log plot in Fig. 3(b). It took over 
a day to obtain one of the relaxation curves in the figure. As lowering the temperature, the 
relaxation curve loses its curvature and tends to become a straight line as can be seen in 



We focus attention on the power-law index obtained by fitting the relaxation curves to the 
function r", where a is a constant. Figure 4 shows the temperature dependence of a, which 
stays constant in the high temperature region, shows a kink at around 1.5 K, and increases 
steeply to unity at low temperatures. This kink temperature is nearly the same as the onset 
temperature of the slush state, where the resistance starts to decrease steeply toward the 
residual low value |l|. The non- unity value in a and its anomalous temperature dependence 



Fig. 3 (a). 
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are characteristic of the present system, whereas the other K-type organic superconductors, 
such as /t-(BEDT-TTF)2Cu[N(CN)2]Br, which has less two-dimensionahty and hence does 
not show any These facts indicate that the resistance 

drop observed by Sasaki et al. corresponds to a microscopic change in the vortex state. 
The sharpness in the kink of the power index a suggests that the change is not a crossover 
but corresponds to a remnant of the vortex glass-hquid phase transition, which has been 
obscured by quantum fluctuation and displaced to lower temperatures 

Generally, power-law-type relaxation curves with non-unity a come from the distribution 
in Ti, and are often seen in disordered systems, such as non-magnetic systems containing 



dilute relaxation centers 



18| . In these systems, each nuclear spins feel the field fluctuation 



with a different amplitude depending on the distance from the nearest neighboring relaxation 
center, and hence has a different Tj~^, which is proportional to the squared fluctuation 
amplitude. Since the present system is in the clean limit and free from magnetic impurities, 
we consider that the relaxation center has its origin in, for example, crystal dislocations 
which has an extremely low spatial concentration. The slight dependence of a in the liquid 
state on the sample jjl and {12 may be rooted in the difference in their concentration in 
each crystal. The vortices in the liquid state contribute little to the nuclear spin relaxation, 
because they are considered to have a field fluctuation spectrum weighted at the frequency 
much higher than the NMR Lamor frequency. Upon entering the slush state from the liquid 
state, there appear fragments of vortex glass in the liquid. These fragments, which waft 
slowly in the liquid is expected to drive the nuclear relaxation. As the concentration of the 
fragments increases with decreasing temperature, all the proton nuclear spins feel the field 
fluctuation with nearly the same amplitude, and will show a simple exponential recovery, 
that corresponds to a=l. 

Finally, we comment on the effect of the nuclear spin diffusion on the relaxation. The 
appearance of power-law-type relaxation curves indicates that the fast nuclear spin diffusion 
is absent in this system. In the presence of fast nuclear spin diffusion, the normal electrons 
in the vortex core must dominate the nuclear spin relaxation, which is expected to be 



proportional to the applied field, as is often observed in high-Tc cuprates|16l. Il7l|. Note that 
in this system, on the contrary, is suppressed in the higher field. We can estimate 
the mean distance between dislocations that act as relaxation centers, from the condition 
that the nuclear spin diffusion does not propagate from one relaxation center to its nearest 
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FIG. 4: Temperature dependence of a, the power law index of the relaxation curve, the time 
evolution of which obeys to the power law r". The dashed curves are guides for the eyes. 

neighboring center within a period of time r ~ 2 x 10^ s in which power-law-type relaxation 
is observed. The lower limit of the mean distance is given as 

i ~ ~ v^(L2/T2)r ~ 1 X 10^ - 3 X lO^A 



(2) 

where L ~ 5A is the mean distance between BEDT-TTF molecules, T2 ~ 30-100 /is, the nu- 
clear spin-spin relaxation rate between the two proton nuclei belonging to adjacent two 
BEDT-molecules. This estimation is consistent with the report by Shubnikov-de Haas 
effect [19], supporting the validity of our interpretation of the power-law-type relaxation. 



IV. SUMMARY 



We have investigated NMR spectra and relaxations in the quasi-two-dimensional super- 
conductor K-(BEDT-TTF)2Cu(NCS)2. The vortex glass-liquid transition under a low field 
of around 0.75 T was confirmed from the diverging T^^ and the spectral broadening. Under 
high field around 4.6 - 5 T, the temperature dependence of the power index a of the nuclear 
spin relaxation curve showed a kink at the vortex liquid-slush transition temperature, below 
which, a approached unity. This behavior of a is understood in terms of an enhancement of a 
short range order of vortices on entering into the slush state, indicating the first microscopic 
observation of the vortex slush. 



8 



Acknowledgments 



This work was supported in part by a Grant-in- Aid for Scientific Research on the Priority 
Area "High Magnetic Field Spin Science in lOOT" from MEXT, and by the Kurata Memorial 
Hitachi Science and Technology Foundation and Saneyoshi Scholarship Foundation. A part 
of this study was performed at the High Field Laboratory for Superconducting Materials, 
Institute for Materials Research, Tohoku University. 



[1] T. Sasaki, T. Fukuda, T. Nishizaki, T. Fujita, N. Yoneyama, and N. Kobayashi, W. Biberacher, 

Phys. Rev. B66, 224513 (2002). 
[2] T. Sasaki, W. Biberacher, K. Neumaier, W. Hehn, and K. Andres, T. Fukase, Phys. Rev. B57, 

10889 (1998). 
[3] R. Ikeda, J. Phys. Soc. Jpn. 72, 2930 (2003). 

[4] T. Nishizaki, T. Sasaki, T. Fukase, and N. Kobayashi, Phys. Rev. B54, R3760 (1996). 
[5] S. M. Dc. Soto, C. P. Slichter, H. H. Wang, U. Geiser, J. M. WilUams, Phys. Rev. Lett. 70, 
2956, (1993). 

[6] H. Mayaffre, P. Wzietek, D. Jerome, S. Brazofskii, Phys. Rev. Lett. 70, 2956, (1993). 
[7] A. Van-Quynh, C. Berthier, H. Mayaffre and P. Segransan, P. Batail, Phys. Rev. B59, 12064, 
(1999). 

[8] S. L. Lee, F. L. Pratt, S. J. Blundell, C. M. Aegcrtcr, P. A. Pattendcn, K. H. Chow, E. M. 

Forgan, T. Sasaki, W. Hayes, and H. Keher, Phys. Rev. Lett. 79, 1563, (1997). 
[9] F. L. Pratt, S. J. Blundell, T. Lancaster, M. L. Brooks, S. L. Lee, N. Toyota, T. Sasaki, 

Synthetic Metals 152 (2005) 417-420. 
[10] A. G. Redfield, Phys. Rev. 162, 368, (1967). 

[11] T. Takahashi, T. Tokiwa, K. Kanoda, H. Urayama, H. Yamochi, G. Saito, Synth. Metals 27, 
A319 (1988). 

[12] T. Takahashi, K. Kanoda, K. Akiba, K. Sakao, M. Watabe, K. Suzuki, G. Saito, Synth. Metals 
42, 2005 (1991). 

[13] T. Takahashi, T. Tokiwa, K. Kanoda, H. Urayama, H. Yamochi, G. Saito, Physica C153-155, 
487, (1988). 



9 



[14] T. Takahashi, K. Kanoda, G. Saito, Physica C185-189, 366, (1991). 

[15] T. Giamarchi, P. Le Donssol, Phys. Rev. B55, 6577 (1997). 

[16] K. Magishi, Y. Kitaoka; K. Asayama, Phys. Rev. B53, 8906 (1996). 

[17] K. Ishida, Y. Kitaoka, K. Asayama, Solid State Comm. 90, 563 (1994). 

[18] M. R. McHenry, B. G. Silbernagel and J. H. Wernick: Phys. Rev. Lett. 27 (1971) 426; Phys. 
Rev. B5 (1972) 2958. 

[19] T. Sasaki, T. Fukuda, N. Yoneyama, and N. Kobayashi, Phys. Rev. B67, 144521 (2003). 



10 



